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Abstract 

In this note we extend integrability conditions for the symmetric stretch tensor U in the polar decom¬ 
position of the deformation gradient V<£ = F = RU to the non-symmetric case. In doing so we recover 
integrability conditions for the first Cosserat deformation tensor. Let F = RU with R : SI C R 3 —> SO(3) 
and U : fi C R 3 —> GL(3). Then 


A := R Grad R = Anti 


|h(Curl U) T - ±tr(U(CurlI/) t )i]eF 


det U 


giving a connection between the first Cosserat deformation tensor U and the second Cosserat tensor Ji. 

(Here, Anti denotes an isomorphism between R 3x3 and 6o(3) := {21 € R 3x3x3 | 21 . M g 50 ( 3 ) \/tt € R 3 }.) The 
formula shows that it is not possible to prescribe U and A independent from each other. 

We also propose a new energy formulation of geometrically nonlinear Cosserat models which completely 
separate the effects of nonsymmetric straining and curvature. For very weak constitutive assumptions (no 
direct boundary condition on rotations, zero Cosserat couple modulus, quadratic curvature energy) we show 
existence of minimizers in Sobolev-spaces. 

Keywords: Cosserat continuum, geometrically nonlinear micropolar elasticity, integrability conditions, com¬ 
patibility conditions, extended continuum mechanics, strain and curvature measures 
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1 Introduction 

Classical continuum mechanics considers material continua as simple point-continua with points having three 
displacement-degrees of freedom, and the response of a material to the displacement of its points is characterized 
by a symmetric Cauchy stress tensor presupposing that the transmission of loads through surface elements is 
uniquely determined by a force vector, neglecting couples. But such a model is insufficient for the description 
of deformable solids with a certain microstructure, such as cellular materials, foam-like structures (e.g., bones), 
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periodic lattices. That is the point where the geometrically nonlinear Cosserat model comes in. 

In order to unify field theories embracing mechanics, optics and electro-dynamics by means of a common prin¬ 
ciple of least action (Euclidean action), the brothers Francois and Eugene Cosserat tried to find the correct 
general form of the energy for a deformable solid (see [5]). In contrast to the classical theory any material 
point in a Cosserat continuum (also called micropolar continuum) is allowed to rotate without stretch, like an 
infinitesimal rigid body. Thus, in addition to the deformation field ip an independent rotation field R is needed 
to fully describe the generalized continuum. The Cosserat model has therefore 6 degrees of freedom (3 for the 
displacements and 3 for the rotations of points), in contrast to the classical elastic continuum, which has only 
3 degrees of freedom. 

From a historical perspective, the kinematical model of such generalized continua was introduced by the Cosserat 
brothers in [5], but they did not provide any constitutive relations. Fifty years later, Ericksen and Truesdell ;8j 
have developed this idea and have drawn anew the attention to the theory of generalized continua. In the 1960s, 
several variants of the theory of media with microstructure were proposed by Toupin [44J , Mindlin [23] , Eringen 
[5] and others. The micropolar (or Cosserat) continua can be viewed as a special case of the nricrostretch con¬ 
tinua and of the micromorphic continua (see Eringen [lU] for the foundation of this theory). In the last decades, 
the generalized theories of continua have proven useful for the treatment of complex mechanical problems, for 
which the classical theory of elasticity is not satisfactory. In this respect, we mention for instance the works of 
Lakes [d [TB] [TS], Neff [29], Neff and Forest [31] . among others. 

The mathematical problem related to the deformation of Cosserat elastic bodies has been investigated in many 
works. In the linearized micropolar elasticity the existence and the properties of solutions have been studied 
by Ie§an .12, 13, 14] and by Neff [27] [15] under some weaker assumptions, among others. For the geometrically 
nonlinear micromorphic model, the first existence theorem based on convexity arguments was presented by Neff 
[25] [2T i ;2BJ. Another existence result for generalized continua with microstructure was proved in [2T] , under 
certain convexity assumptions. Adapting the methods from mm, Tambaca and Velcic m have derived an 
existence theorem for nonlinear micropolar elasticity, for a general constitutive behaviour. 

In our paper we consider a new energy formulation of geometrically nonlinear Cosserat models which completely 
separates the effects of nonsymmetric straining and curvature. 

In Section[2]wc introduce the strain measures and curvature strain measures associated to the classical nonlinear 
Cosserat model. The next section presents the expression of the elastically stored energy density and recalls 
the main assumptions on the constitutive coefficients. 

In Section [3] we provide an easy and self contained approach towards compatibility conditions in micropolar 
media. Starting point is Vallee’s formula (Theorem[2]) providing compatibility relations generated by the classical 
polar decomposition. As it turns out, this formula can be generalized to the situation F = R-U, where R is still 
orthogonal, but U need not be positive definite and symmetric. Our final compatibility result is already known. 
Indeed, Teresi and Tiero (cf. [43] , eq. (22)) use a much more abstract setting with an extremely short algebraic 
proof as compared to the classical, componentwise approach (cf. [15]). Their result seems not to be very well 
known and it has not been immediately connected to the classical compatibility result for the symmetric stretch 
tensor. In our presentation we choose a compact direct tensor notation which is a middle between the two 
extreme cases above. 

The advantage of Teresi and Tiero’s analysis is, however, to make the structure more transparent and to high¬ 
light the appearance of the curl of the first Cosserat deformation tensor. This dependence is not easily grasped 
in the compatibility conditions given in the classical approach (see [If)] and Eringen’s book [lU], eq. 1.7.4). 

Having proven these relations between dislocation density K, strain measure U and its derivative in the form 
of CurlC/, in Section [5] we finally give the new energy formulation and briefly outline a proof for the existence 
of minimizers. 
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1.1 Notation 


In this paper, we denote by R 3x3 the set of real 3x3 second order tensors, written with capital letters 
and the Einstein summation convention over repeated indices is used. For a, b £ R 3 we let (a, b) denote the 
scalar product on R 3 with associated vector norm ||a || 2 = (a, a). Similarly, the standard Euclidean scalar 
product on R 3x3 is given by (X, Y) = tr(XY T ), and thus the Frobenius tensor norm is ||X || 2 = (X, X). The 
identity tensor on R 3x3 will be denoted by 1, so that trX = (X, 1). We adopt the usual abbreviations of 
Lie-algebra theory, i.e., so(3) := {X £ R 3x3 \X T = —X} is the Lie-algebra of skew symmetric tensors and 
sl(3) := {X G R 3x3 | tr(X) = 0} is the Lie-algebra of traceless tensors. Sym(3) = {x £ R 3x3 | X = X T } denotes 
the set of symmetric matrices. For all vectors £,77 £ R 3 we have the tensor product (£ © r])ij = rjj and e is 
the Levi-Civita tensor, also called the permutation tensor or third order alternator tensor, given by 



otherwise. 


is an even permutation of (1,2,3) 
is an odd permutation of (1,2,3) 


( 1 ) 


For all X £ R 3x3 we set sym X = 4(X T + X) £ Sym(3), skewX = 4(X — X T ) £ so(3) and the deviatoric part 
devX = X — i(trX)-l £ s((3) and we have the orthogonal Cartan-decomposition of the Lie-algebra gl(3) 


gl(3) = { sl(3) fl Sym(3) } © so(3) © R-1, 
X = dev symX + skewX + — (tr X) ■ 1. 


Here, for 


A = 


we consider the operators axl: so(3) 



f 0 

-0.3 

a 2 \ 

1 ( 

( a i\ 

1 

axl 

03 
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0 ) 
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l «3 / 



( 2 ) 


(3) 


Aij — e?yA:(axl .' 1 )/,. —. (anti (axl /l)), 


ij, 


for all v £ . 


(axlH)^ — 2 f ijk 


(4) 


Furthermore the axial operator Axl : ©o(3) 
defined by 


x3 associated to an so(3)-valued third order tensor 21 is 

(5) 


where 6o(3) := { 21 € 


(Axl 21.u) x v = (21 .u)v, u,v £ I 3 , 

3x3x3 | 21 .u £ so(3)} and we define Anti : R 3x3 —> ©o(3) to be 
&ijk = ~ e iji (Axl &)ik =■ (Anti (Axl©)) , jk . 

Then, for A £ R 3x3 we have 

e . Anti A — (e^fe c* © ej © e&) . ( e rs tA^ u e r © e s © e u ^ — eij^e rs tA^ u Sjrdks © e 
— eijk£jktAt u et © e u 2 SnAt u © e u — 2 A$ u © e u — 2 A 
and thus, applying this to A = Axl 21, 

(Axl 21)= ~(c ■ 21 )ij = —ei rs 2l rS j . 


( 6 ) 


(7) 

Here, the double dot product “ : ” of two third order tensors A = Aijk e* © e,- © ek and B = Bijk e, © e 7 - © ek 
is defined as A : B = Ai rs B rs j © ej . The cross product “x” of two second order tensors is calculated with 
help of the rule (a © b) x (c © d) = a © (6 x c) © d , which holds for any vectors a, b , c and d. 

Throughout this paper (when we do not specify else) Latin subscripts take the values 1, 2, 3. Typical conventions 
for differential operations are implied such as comma followed by a subscript to denote the partial derivative 
with respect to the corresponding cartesian coordinate (e.g., f,i = ). 
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2 Strain and curvature measures for the classical Cosserat model 

We consider an elastic body B = {Vk} as a set of particles Vk, which in the reference configuration occupies a 
bounded domain C R 3 with Lipschitz boundary dfl. Let Ox 1 X 2 X 3 be a Cartesian coordinate frame in R 3 and 
ei the unit vectors along the coordinate axes Oxi . Then any material point Ai in Vk may be identified by the 
spatial location p(x) of the centroid x = (£ 1 , 22 , £ 3 ) of Vk and its relative position vector £( 2 ). Furthermore 
the behaviour of the body under the action of some external loads can be described by the deformation field 

ip: —> R 3 , p(x) = p(x 1 , 22 , 23 ) 

and a second order tensor field R (also called microrotation tensor ) according to 

R: O —>-SO(3), £( 2 ) = R{xi, x 2 , 2 3 )£ 0 


with some fixed £0 £ K 3 - 

We identify the second order tensor R with the 3x3 matrix of its components R £ R 3x3 and the vector p with 
the column-vector of its components p £ R 3 . Notice that R is a proper orthogonal tensor, i.e. R £ SO(3). The 
three columns of the matrix R will be denoted by d\ , d 2 , d 3 . They are usually called the directors and can be 
interpreted as an orthonormal triad of vectors {di} rigidly attached to each material point, describing thus the 
microrotations. Using either the direct tensor notation or the matrix notation, we can write respectively 

R = di®ei or R = (di | d 2 \ d 3 ) 3x3 . 

In the reference configuration tt the directors attached to every material point ( 21 , 22 , 23 ) are taken to be { } 

and after deformation they become {^( 21 , 22 , 23 )}. If we denote by po(xi,x 2 ,X 3 ) := Xi e* the position vector 
of points in the reference configuration f 2 , then the displacement vector field is u = p — po . 

To introduce nonlinear strain measures let us now consider a material point at position 2 £ Vk and let 2 + h 
be a neighboring material point within Vk- Then we get 

p(x + h) — p(x) = Grad p{x) h + o(h), 

which provides us with a primary measure of deformation, the so called deformation gradient 

F = S/p = gradtp = p,i®ei = (p, 1 | p , 2 \ P,3 ) 3x3 - ( 8 ) 

With this the Cosserat-brothers discussed an energy functional of the form 

I(p,R)= f W(F, R, Grad R) dV — II( R), where 

Jn (9) 

II( p, R) = the external loading potential 

and postulated invariance under Euclidean transformations (today also known as the principle of frame- 
invariance), which states that the energy density has to be left invariant under rigid rotations 

W{QF, QR, QGradf?) = W{F, R, Grad R) for all Q £ SO(3). 

_ fj 1 

With Q = R we get the reduced representation (cf. [SJ[7]) 

W(F, f?,Gradf?) = W(R T F, R T R, R T Grad R) = W(lf F,if Gr&AR ), 

which motivates to consider the first Cosserat deformation tensor (also called the nonsymmetric right 
stretch tensor, see the original Cosserat book[5], p. 123, eq. (43)) 

U = if F = (a 0 di)(p,j 0 ej) = {di, p,j ) e, ® ej . (10) 


4 


and the second Cosserat deformation tensor (also called the third order curvature tensor , see the original 
book [5], P- 123, eq. (44)) 


.ft = i? T Gradf? = R T (R, k 0 e k ) = (R T R, k ) 0 
— ( di , dj^k ) ei 0 €j 0 ek — (R R,i \ R' R ,2 

as proper deformation measures. Although A is a third order tensor, 
components, since R R, k is skew-symmetric (k = 1,2,3). 

Remark 1 . We write the above tensors R and U with superposed bars in order to distinguish them from the 

factors R and U of the classical polar decomposition F = RU, in which R is orthogonal and U is positive 

definite, symmetric and which is a standard notation in elasticity. 

Then, according to the classical Cosserat theory, we define the relative Lagrangian strain measure for stretch 

by 

I = U - 1 3 = (( di , ip,j ) - Sij ) ei <g> ej , (12) 

where Sij is the Kronecker symbol and I 3 := 0 e* is the unit tensor in the 3-space and for the Lagrangian 

strain measure for orientation change (curvature) we use the second Cosserat deformation tensor. 

If we take the transpose in the last two components of A, then we obtain the third order curvature tensor A 
which was used in [23 [23 [ 11 ] ; 

~ 2.3 

— (dt, dk,j ) e^ 0 ej 0 e k — (e^ 0 dj) (dk.j 0 fo) 0 e k (13) 

= i? T (Graddfc) 0 et = (R T Grad g?i | i? T Gradc ?2 | R T Grad g ?3 ) . 

In order to avoid working with a third order tensor for the curvature, one can replace A by a second order 
curvature tensor. This can be done in several ways. Indeed, a first way is to consider the axial vector of the 
skew-symmetric factor R R,k in the definition m of A and to introduce thus the second order tensor 

T = axl(i? T i?, fe ) 0 e fc . (14) 

The tensor T is a Lagrangian measure for curvature and it is frequently called the wryness tensor in the literature 
(see e.g., [22]). The relation between T and M. can be expressed with the help of the third order alternator tensor 

e = - 1 3 x 1 3 = djk ei 0 ej 0 e k 

in the form (since axl S = — 5 e : S, for any skew-symmetric second order tensor S) 

T = — — e : A and A = I 3 x T = — e T. (15) 

As an alternative to the wryness tensor T given by (fTTl) , one can also use the Curl operator instead of Grad in 
the definition (HD of A. In doing so, we define the dislocation density tensor K by 

K = Tf Curl R = -R T (R,i x a) . (16) 

The dislocation density tensor K presents some advantages for micropolar and micromorphic media, as it can 
be seen in El [21 [31 US] ■ Here, the Curl operator for vector fields v — Vi ei has the well-known expression 

curl v = tijk Vjj e k = —v,i xe* , 

while the Curl operator for tensor fields P = Pjj ei 0 e 3 is defined as 

Curl P — ej rs Pis,r ei 0 ej — P,i x ei , (1^) 


e k 

I R t R ,3 ) , 


( 11 ) 


it has, in fact, only 9 independent 
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In other words, Curl is defined row wise as in mm]: the rows of the 3x3 matrix Curl P are the 3 vectors 
curl Pi (i = 1,2,3), respectively. Note that sometimes other definitions of Curl are used, for example in 1 451 . 
The Curl operator there, which we will denote by Curl, is defined by the requirement (CurlP)u = curl(Pu), 
entailing that ChM(P T ) = (CurlP) T 

We describe next the close relationship between the wryness tensor T, the dislocation density tensor K and the 
second Cosserat deformation tensor. One can prove by a straightforward calculation that the following relations 
hold m 

-K = r T -(trr)l 3 and - T = R T - ^ (tiK) 1 3 , (18) 

as with USD 

£ = -eT = eK T -^(trK)e. (19) 

For infinitesimal strains this formula is well-known under the name Nye’s formula , see [34] . and there (—T) is 
also called Nye’s curvature tensor (cf. [37])- From (fT51) we deduce the following relations between the traces, 
symmetric parts and skew-symmetric parts of these two tensors 

tr K = 2trT, skew AT = skewT, devsymA' = — devsymT. (20) 

In view of m GOD, we see that one can work either with the wryness tensor F or, alternatively, with the 
dislocation density tensor K , since they are in a simple one-to-one relation. 

3 Integrability conditions for the symmetric stretch tensor U 

In classical continuum mechanics, the deformation of a continuous medium can be characterized by the Green- 
Lagrange strain tensor, which is defined by E = \ (P T F — 1 3 ) = \ (\7ip T V(p — 1 3 ) . Conversely, if a positive 
definite symmetric matrix C = I 3 + 2E is given in terms of x, it would be interesting to know under what 
conditions the system 

XLp T \7<p = C (21) 

of first order partial differential equations is liable to be solved. In mm® the authors have shown proper 
compatibility conditions in the case of a classical continuum undergoing large deformations, which yields the 
compatibility relation 

(Curl A) t + Adj ^4 = 0, (22) 


where Adj A is the adjugate of A and with 

A = ( ([/(Curl U) T - \ tr([/(Curl U) T ) 1 3 ) Uj . 

In addition (cf. [35]), C. Vallee gets the following useful result for the symmetric stretch tensor U = C 1 / 2 . 

Theorem 2. Let LI C R 3 be a domain. Suppose ip: 12 —► R 3 is an twice continously differentiable deformation 
field of an elastic body with polar decomposition 

S7<p{x) = R(x) U{x), (23) 

where R: Ll —> and U: Ll —>■ Psym(3) is positive definite, symmetric and differentiable. Then for any x £ Q 
and u, v £ R 3 we have 

[(P T GradP)(x).u] .v = ^ t ^ ^ (|p(a:)(CurlP(x)) T — ^ tr(t/(x)(Curl [/(a;)) T )l| U(x)iij x v (24) 


6 




Alternatively to (liHl) . we have for x £ Ll 

Axl (R t Grad R) (x) = —^—{[/(x)(Curl[f(x)) T - Lr([/(x)(CurlC(x)) T )ll C(x). (25) 

det U(x) l 2 ) 

Remark 3. Let C R 3 be a domain. Theorem H is also well-known for the case of infinitesimal deformations 
ip: f l —> R 3 (see m in which the polar decomposition V ip = F = RU turns into 

1 + Vti = (1 + skewV u + .. .)(1 + sym Vu + ...) (26) 

in the neighbourhood o/l, where skewVu is the infinitesimal continuum relation and sym Vu is the infinitesimal 
stretch. Thus 

Vu = sym Vu + skewVu (27) 

is the linearized polar decomposition and it holds 

V axl(skewVu) = [Curlsym Vu] T — — tr ^( Curlsym Vu) 2 )l = [ Curlsym Vu] T , (28) 

which is the linearization analogue of (1251) . 


4 Integrability conditions for the non-symmetric first Cosserat de¬ 
formation tensor U 


In [IB] the authors obtained a compatibility condition for the right Cauchy-Green tensor, using homology and 
homotopy group techniques and in [38j [40] one can find compatibility conditions for rotations. In a more direct 
approach we get the following result similar to Theorem [5] 

Theorem 4. Let Ll C R 3 be a domain and < p: f l — > R 3 be a twice continuously differentiable function. Suppose, 
for any x € fi, 

V<p(x) = R(x) U(x), (29) 

where R: f l —> 0 ( 3 ) and U: fI —> GL(3) are continuously differentiable. Then, for any x € LI 

Axl ( R T Grad R) (x) = -^—|l7(x)(CurU7(a;)) T — - tr(?7(x)(Curl ?7(o;)) T )l} U(x). (30) 

det U(x) 1 2 J 

Remark 5. Similar to Remark [5| we obtain the formula 

V axl (A) = ^ Curl(Vu — A)j — i tr f ^ Curl(Vu — A)^ j ■ 1 


as linearization of E3D, where Vu — A is the linearization of the first Cosserat deformation tensor and V axl(A) 
is the linearization of the second Cosserat curvature tensor. 

As partial converse, we obtain formula (1221) as condition for the existence of a deformation tp admitting the 
decomposition JUD also in the case of not necessarily symmetric tensors U: 

Theorem 6 . Let LI C R 3 be a simply connected domain. Let U: LI —> R 3x3 be continuously differentiable with 
values in GL(3). Then there exists a deformation < p: —> R 3 and R: LI —> 0(3) satisfying 


V</?( x) = R(x)U(x) for all x £ LI, 

provided that 

A ■- -2= [ F(CurlI7) r - l tr(F(CurlI7) T )l ) U (31) 

det U V 2 J 

is continuously differentiable and satisfies the compatibility condition 

(Curl A(x)) T + Adj A(x) = 0 for every x £ LI. 

Furthermore this deformation is unique up to a rigid body deformation, that is multiplication with a constant 
rotation and addition of a constant vector. 
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The representation given by Theorem [5] makes it possible to express the Cosserat tensors in terms of U and 
Curl U in the following way: 

Consequence 1. Let Q C R 3 be a domain and ip: Ll —> R 3 be a twice continuously differentiable function. 
Suppose, we have the decomposition 


SRp(x) = R{x) U(x) for any x £ LI, 


(32) 


where R : f l —> 0(3) and U : f l —> GL( 3) are differentiable. 


Then, we have 


A = R Grad R = Anti 


K = —A : e = — Anti 


U(x)(Cm\U(x)) T — \ tr(J 7 (a;)(Curl[/(x)) T )l 3 U(x) 


det U 


U(x)(Cur\U(x)) T — ^tr({ 7 (x)(Curl[/(x)) T )l 3 U(x) 


det U 


: e. 


(33) 

(34) 


Proof. Equation (1551) directly follows from the definition © of Anti and (1301) . Formula (1551) hence is a conse¬ 
quence of the identity K = —A : e, which has been proven in [30|- □ 

Another consequence of Theorem [I] is that the Curl of U can be written as function of U and K as follows. 
Note, however, that it is to be expected that the conditions concerning regularity and structure (orthogonality 
for R and gradient-structure of R U simultaneously) make extending this consequence to less regular functions 
more difficult than the usual trivial approximation by sequences of mollified functions. 


Consequence 2. Let O C R 3 be a domain and <p: O —> R 3 be a twice continuously differentiable function. 
Suppose we have the decomposition 


V</?(x) = R(x) U(x) for any x £ LI , 

where R : f l —> 0(3) and U : f2 —> GL( 3) are continuosly differentiable. Then, we have 
(Curl U) T = Adj U (TU- 1 — tr (vU J.3^ , that is 


(Curl U) T = Adj U 
Proof. Again we define: 


tr(AT) (U ~ 1 - tr (U _1 ) Is) - (W _1 - tr (a^C/ ) 1 3 ) 


A := -= 

det U 


1 r _ 1 


[/(Curl U) 1 - -tr([/(Curl[/) i ) 1 3 


U, 


(35) 

(36) 

(37) 


as in (ED, so that 0 readily shows e : A = — 2AxlA = — 2 A according to ([30]) . Thus we have (using ([T5)l l 
T = — : A. = A and (1551) follows from Lemma ITT! below or, more precisely, from the first step of the proof of 

Theorem QJ that is (H51) . Therefore also (1571) . being the combination of (1551) and (1151) . is proven. □ 


4.1 Proof of the Theorem 

For the sake of notational convenience let us introduce the following notation in order to interpret the determi¬ 
nant of a matrix as function of its column vectors: 

Definition 7. In the following, let vol: R 3 x R 3 x R 3 —> R denote the 3-linear alternating 3-form on R 3 defined 
by 

vol(it, v, w) = det A, u,v,w GS 3 

where A = (u|u|w). 













A direct consequence of the definition then is, that for A £ R 3x3 and u,v,w £ M 3 we have (cf. Appendix. 
Lemma E3 


vol(Au, Av, Aw) = det A vol(w, v, w) and vol(Au, Av, w) = vol(«, v, (Adj A)w), (38) 

where Adj A denotes the adjugate of A. 

If the matrix is applied to only one argument of vol, we obtain the following formula (cf. Appendix, Lemma 

EH): 


vol(Aw, v, w) + vol(u, Av, w) + vol(w, v, Aw) = (tr A) vol(it, v, w ) (39) 

for all u, v, w £ R 3 and A £ R 3x3 . 

Lemma 8 . Let A,B£ R 3x3 , det B ^ 0, let h,k £ R 3 . Then 

(Ah) x ( Bk) - (Ak) x ( Bh) = [tr(AB -1 )(Adj B) T - ((Adj B)AB~ 1 ) T ] h x k. 

Proof. Let us consider the scalar product of the left-hand-side with arbitrary »eR 3 . Since (axb,v) = vol(a, b, v) 
for all a, b, v £ R 3 (cf. Appendix, Lemma fliTl) we have 


L := (Ah x Bk — Ak x Bh, v) = vol (Ah, Bk, v) — vol(Afc, Bh, v) 

= vol (BB~ X Ah, Bk, v) — vol (BB~ X Ak, Bh, v), 


which, by (138L equals 

vol (B~ 1 Ah, k, Adj Bv) — vo\(B~ 1 Ak, h, Adj Bv) 

and after switching two arguments of the second vol and adding and substracting 
vol(/i, k, B~ l A Adj Bv), we infer 

L = vol (B~ 1 Ah, k, Adj Bv) + vol(/i, B~ 1 Ak, Adj Bv) 

+ vol (h, k, B~ X A Adj Bv) — vol (h, k, B~ 1 A Adj Bv). 


Here (1551) can be applied to combine the first three terms and since (a x b, v) = vol (a,b,v) (see (|55|) L we get 

L = tr(B~ 1 A) vol (h, k, Adj Bv) — vol(h, k, B~ X A Adj Bv) 

= tr (B^ 1 A)(h x k , Adj Bv) — (h x k, B~ x A Adj Bv) 

= tr(H _ 1 A)(Adj B T h x k,v) — ((5 -1 A Adj B) T h x k,v) 

= ([tr(H _ 1 A) Adj B t - (B- 1 A Adj B) T ] hxk,v). 

Due to the arbitrariness of v, this proves the claim if we take into account that tr (B~ 1 A) = tr (AB~ l ). □ 

Lemma 9. Let A £ R 3x3 . Then, for any h,k,v £ R 3 , 

((Adj A) T (h x k)) x v = Ah x (Ak x v) — (Ak x Ah) x v. 

Proof. Let us recall that, for any a,b,c £ R 3 , a x (b x c) = (a, c)b — (a, b)c. Using this rule twice, for arbitrary 
h,k,v £ R 3 we obtain 


Ah x (Ak x v) — Ak x (Ah x v) = (Ah, v)Ak — (Ah, Ak)v — (Ak, v)Ah -I- (Ak, Ah)v 

=(Ah, v)Ak — (Ak, v)Ah = v x (Ak x Ah) = (Ah x Ak) x v = ((Adj A) T (h x k)) x v, 

where the last equality results from the fact that for any h, k £ R 3 and A £ R 3x3 , Ah x Ak = (Adj A) T (h x k) 
due to the second part of dSHD- □ 

Lemma 10. Let U : —> R 3x3 be differentiable. Then for any x £ LI, h,k £ R 3 , we have 

(Grad U(x).h)k - (Grad U(x).k)h = Cm\U(x).(h x k). 
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Proof. According to the Definition we have 

Curl U (%) .{h X k'j — ( t /,2 (x) ^ {fir e x X k s e s ) — {fijkGmk,j (x) ® &i) ffrsthskt e r ) 

— £ijk£-rstUmk,j {x^)h s ktSi r G"rn — ^ijk^-istU rn k,j{x > )h s kt e m 

= £jik£istUmk.j{x)h s kt e m = {Sj s 5kt djt5ks)U rn k,j{x)h s kt e m 

= ( Umt,s{U ms ^{x^h s kt e m 

and on the other hand 

(Gradf/(x).h)fc — ( Grad U{x).k)h = (( U,k (x) Cg> ek)hi ei)k r e r — (( U,k (x) Cg> ek)k r e r )hi e; 

— ( J-^ij,k{x)hlSlk Cj ® Ej'jkr € r {Uij^k{x')k r $kr ® El 

— {Uij,l(x)hik r 5j r Uij^ r {x')k r hiSji'^ Ei 

— (Uir,l{x) Un^ r {xf^hlk r El 

and therefore the lemma is proved. □ 

Lemma 11. Lot R: S2 — > R 3x3 bs a diffErsntiable function with values in 0(3). Then for any x £ fl, h £ R 3 , 

A.h = R T {x) Grad R{x).h £ R 3x3 is skew-symmetric (40) 

and there is a map A: SI x R 3 —> R 3 , linear with respect to its second argument, such that for any a; £ SI, 
h £ R 3 , v £ R 3 , 

(Grad R{x).h)v = R(x)A{x, h) x v . (41) 

Proof. Because R{x) £ 0(3) for any x £ SI, 

I 3 = R T (x)R{x) 

holds true for any x £ LI and differentiating leads to 

0= [Gradl?(x)./i] T .R(x) + R{x) T GradR{x).h = 2sym{R{x) T GradR{x).h) 

for arbitrary x £ SI, h £ R 3 , which shows that the symmetric part of R{x) T Grad R{x).h vanishes and implies 
pnii . By definition 03> R (x) Grad R{x).h therefore can be expressed as 

(R {x) Grad R{x).h)v = A{x, h) x v, 

where A{x,h) depends linearly on h. Now, equation (14111 follows directly upon multiplication by R{x). □ 

In order to assert solvability of differential equations, let us recall the Theorem of Frobenius (see 071 Thm. 4.E, 
p. 167]): 

Lemma 12. Let X,Y be Banach spaces, U C X xY open and F: U —> L{X,Y) a continuously differentiable 
function. Then the following are equivalent: 

i) For any (xo,yo) £ U, the differential equation Dy{x ) = F{x,y{x)), y{x 0 ) = yo has exactly one C 1 -solution in 
a neighbourhood of{xo,yo). 

ii) D x F{x , y)hk + D y F{ x, y)F{x, y)hk = D x F{x , y)kh + D y F{x, y)F{x , y)kh for all (x, y) £ U and all h,k £ X. 
We immediately take this lemma to its use in the following: 

Lemma 13. Let SI C R 3 be simply connected. Given a differentiable function A: SI — > R 3x3 such that 

(Adj A{x)) T + Curl A{x) = 0 for all x £ S2, 

there exists a continuously differentiable function R: SI —» R 3x3 with values in 0(3) such that 

[Grad R{x).h\v = R{x){{A{x)h) x v) for all x £ 0, h,v £ R 3 . (42) 

This function is unique up to the multiplication by a constant Q £ 0(3). 
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Proof. We want to apply Frobenius’ theorem (Lemma 1T21) with F: 12 x R 3x3 —> L(R 3 ,R 3x3 ) defined by 

[F(x, Y)h\v = Y[A(x)h}v, iefi.be M 3x3 , h, v £ R 3 . 

Thus, the partial derivatives of F are given by 

[D x F(x, Y)hk\v =Y[D(A(x)k)h] x v, 

[[DyF{x,R)Z]S]v=Z{A(x)S)v, 

[DyF(x,Y)F(x,Y)hk\v =F(x,Y)h(A(x)k) x v = Y(A(x)h) x ((A(x)k) x v). 

Condition 2 from Frobenius’ theorem (Lemma 1121) therefore is 

Ah x ( Ak x v) — Ak x (Ah x«) + (Curl A(h x k)) x v = 0 V7i, fc, u, (43) 

which we can obtain as follows: 

(Adj A) T + Curl A = 0 

implies that 

(((Adj A) t + Curl A)(h x k)) xv = 0 for all v £ R 3 , h,k £ R 3 . 

By Lemma [9l this entails 


Ah x (Ak xu)— Ak x (Ah x v) + (Curl A(h x k)) x v = 0 V7i, k, v, 

which is mg) . Therefore, since 12 is simply connected, choosing an arbitrary point Xq £ 12 and yo £ 0(3 ) 
yields a unique solution R to (13211 satisfying R(x o) = yo- Since R T DR(x)h is skew-symmetric for any x £ 
12 and any h £ R 3 by (1321) (it can be expressed as vector product), R 1 R is constant within Q and hence 
R T (x)R(x) = y^yo = 1, so that actually R(x) £ 0(3) for all x £ H. Given a second solution R: —> 0(3), by 
the aforementioned uniqueness property, R(xq)R(xq)~ 1 R = R. □ 

Proof of the theorems. The theorems can be seen as partially converse statements. Let us begin their proofs 
with the part they have in common: 

Lemma 14. Let U : Ll — >• R 3x3 be differentiable and U(x) invertible for every x £ LI and let R: 12 —>• R 3x3 be 
a differentiable function with values in 0(3). Then 

Grad(RTl)(x)(k 1 h) = Grad(f? U)(x)(h, k) (44) 


if and only if 

detU(x)[A(x)U 1 (x) — tr(A(x)U 1 (x))la] = t/(a;)(CurU7(x)) T for any x £ 12, (45) 

where A is defined as in Lemma CUD- 
Proof. The product rule implies that 

Grad (RU)(x)(h, k ) = [Grad R(x).k]U(x)h + R(x)[Gr&dU(x).k\h for any x £ 12, h, k £ R 3 
and calculating the difference between two such expressions with different order of the arguments, we obtain 
Grad(.R U)(x)(h,k) — Grad(i? U)(x)(k,h) 

=(Giad R(x).k)U(x)h + R(x)[GradU (x) .k\h — [Grad R(x).h]U(x)k — R(x) [Gr ad U(x).h\k. 
An application of Lemma [TT] shows that (1331) is equivalent to 

0 = R(x)A(x, k) x (U(x)h) + R(x)[GvadU(x).k\h — R(x)A(x,h) x (U(x)k) — R(x)[GradU(x).h]k 

_ fp 

for any x £ 12, h, k £ R, where A(x, h) is defined as in (1311) . Multiplying this equation by R (x) and re-ordering 
the terms, we equivalently arrive at 

A(x)h x U(x)k — A(x)k x U(x)h = — ([Grad U(x)h]k + [Grad U(x)k\h) 
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for any a; G f2, h, k G I 3 , where we can apply the equalities from Lemmata l8l and flOl to the left-hand-side and 
right-hand-side respectively. This transforms the equation into 

[tr (^A(x)U *(x)^ [Adj U{x)) T — ((Adj U(x))A(x)U 1 (x)) T ]h x k = — Curl U(x)(h x k) 
for all x G T 1, h,k £ R 3 and therefore, by transposition and because h, k were arbitrary, 

Adj U(x)[A(x)U 1 (x) — ti(A(x)U ^x))!^] = (Curl t/(x)) T 
for any x£(l. Multiplication by the invertible matrix U(x) from the left gives 

det U{x)[A{x)U 1 (x) — tr(A(x)(7 1 (x))ls] = 17(x)(Curl{7(x)) T for any x GO, 
which finally asserts that da and dlS]) are equivalent. □ 

We continue with the proof of Theorem 2) 


(46) 


Proof of Theorem As V<p = RU, by Schwarz’ theorem and Lemma [HI we obtain 0 = D 2 ip(x)hk — D 2 ip(x)kh 
and hence 

det U(x)[A(x)U 1 (x)—tr (A(x)U 1 (x))l 3 ] = U(x)(CurH7(x)) r for any x G O. 

Taking the trace on both sides, we infer 

det (7(x)[tr(A(x)f7 *(x)) — tr(A(x)U *(x)) • 3] = tr((7(x)(Cmi[/(x)) T ) 
for any x G O and hence 

1 


tr ^A{x)U *(x)^ 
Inserting this into (l45l) shows 
det U(x)[A(x)U 1 (x) 
which entails 


2 det U (x) 


■ tr(C7 (x)(Curl U (x)) J 


for all x G O. 


2 det U (x) 


tr([/(x)(Curl C/(x)) T )l 3 ] = [/(x)(Curl U{x)) T for x G O, 


-=—tr[[/(x)(Curl U(x)) T ]U(x) for all x G O. 


A{x) = -=—-C/(x)(CurlC/(x)) T C/(x)-=— 

detf7(x) 2 det U(x) 

By definition ja of A finally this translates into 

[Gradl?(x).h]x = R(x) [( ^ 1 , U (x)(Curi U{x)) T U (x) - tr[0(x)(CurlI/(x)) T ]0(x))fc] x v 

for any x G O, h, v G R 3 . 


□ 


Proof of Theorem [6J Since A satisfies the compatibility condition, Lemma fl3l enables us to find R : O — > 0(3) 
such that 

[Gradl?(x).h]x = R(x)((A(x)h) x v) for all xGf2, /i,uGt 3 . 

Furthermore, R is determined uniquely up to multiplication with a constant rotation. 

The definition of A ensures that 


AU 1 = — (U (Curl U) T - itr(C/(CurlC/) T )l) 
det U 2 


and hence 




det U (AU -tr (AU )1) =C/(CurH7) J - - tr([/(Curl U) 1 )1 - tr 


U (Curl U) T - -tr([/(Curl[/) T )l 


1 


=U (Curl U) T - i tr(t/(Curl U) T ) 1 - tr(C/(Cuii U) T ) 1 + | tr(I/(Curl U) T ) 1 
=[7(Curl U) T 

that is, condition (1451) from Lemma El is met, so that condition 2 of Lemma El is satisfied for F(x,Y) = 
R{x)U{x) and hence upon choosing xo G fl and p{x o) G R 3 there is a unique mapping p: —> R 3 such that 
\7p = RU. □ 
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5 Existence theorem for a new Cosserat formulation 


We have seen that in fact K can entirely be obtained as a complicated nonlinear function of the strain measure 
U and its derivative in the form Curl U alone; therefore, any formulation directly written as a simple function 
of (U,K) is a function of (U, CurH7), in which, however, the effects of strain U , and generalized curvature 
CurH7 are not separated. Next, we would like to present an energetic formulation which has the possibility 
to additively separate the pure straining part sym (U — 1) and Curl?/ for certain parameter ranges. Such a 
formulation has not been considered previously. We set 


I((p,R)= / W(U,K)dV — H(p,U), where 

Jn 


(47) 

W(tp, R) —W(U, K) = p sym(U — 1) 2 + p c skew(I7 — 

1 ) 2 


+ 9 (det U 2 2) + pL 

||Curl(?7 — 1)|| + pL c ||/\'|| , 

(48) 

2 det U 

n(p,u)= I fpdv 

Jn 




for a given body force 0 ^ € L 2 (fl). Taking into account (1381) and (1371) it is clear that, for sufficiently smooth 

functions, ||Curl([7 — 1)|| = ||CurlC/1| is indeed a complicated, nonlinear function of U and K. 

In (1481) . p c = 0 means that only the symmetric part of the straining is taken into account and L c is an additional 
length scale which makes (|T7|) well-posed. It would be the final goal to show that flUD with W(tp, R ) given by 
(1481) admits minimizers for the limit parameter range p c = 0 and L c = 0 in a suitable Soboloev-space dictated 
by the coercivity of the problem. Unfortunately, currently we cannot dispose of the contribution of ||A~|| by 
setting L c = 0 and thereby completing the separation of the effects of U and those of Curl U. Let us summarize 
what can be obtained for the minimization problem posed by (1471) : 

Case (I): p > 0, A > 0, L c > 0; p c > 0, L c > 0, q > 2 

existence of minimizers along the lines presented in m- 
Case (II): p > 0, A > 0, L c > 0; p c = 0, L c > 0, q > 3 

^ existence of minimizers along the lines presented in pTTl using [58;. 

Case (HI): p > 0, A > 0, L c > 0; p c = 0, L c > 0, q > 2 
^ outlined subsequently 

Case (IV): p > 0, A > 0, L c > 0; p c = 0, L c = 0 
^ yet open 

We assume fixed boundary conditions of place at To C dQ, i.e. <p(x) | r = x. Consistent with this requirement 
we impose zero tangential straining at which means (taking the cross product of a matrix with a vector row 
wise) 

(U — 1) x v = 0, where v denotes the outer unit normal 
4=> (U — l).r = 0 for all tangential directions r at ]?£>. 

The latter condition translates into 

R T \7<p.T = r at To (49) 

and due to p(x) | r = x we have \7p(x).r = r. Thus we require in fact 

R - t \v d = t R \r D =1 (50) 

Since the body force / is non-vanishing, the solution (ip, R) will be non-trivial. 
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5.1 Outline for Case (III) 

We use the direct methods of the calculus of variations and only indicate those arguments which are needed in 
addition to the proof in [SO] for case (I). At first, let us recall the incompatible Korn’s inequality (cf. [2; 135]): 


IlsymPf 


|CurlP| 2 dV > c + (fi) • / ||P|| 2 dV 


(51) 


Applying (l5!l) to P = U — 1, we obtain a uniform P 2 -bound on U — 1 (while only controlling the symmetric 
part locally). Considering infimizing sequences ipk, Rk we obtain (due to the fact that U — 1 is uniformly 
P 2 -bounded), similar to [30], eq. (42): 


h 1 ( n) 

Pk -^ P 


L 2 (Q) 

and ip k - > p, 


Curl U k g(Curi) - Curl U and R T k Vp k = U k U. 

The independent control of Rk is provided by that of ||A'|| 4 , which gives 

Rk R and CurlP fc CurlP. 

On passing to a subsequence, not relabelled, we obtain 



l 2 ( n) 


Pk 

- > 

P, 

Vp k 

l 2 ( n) 

Vip, 

Uk 

L 2 (Cl) 

U = R T Vp, 

CurlPfc 

l 2 ( n)^ 

Curl U, 

Curl Rk 

l 2 ( n) 

Curl R. 


(52) 


(53) 


(54) 


Moreover, 

(det U 2 + (det Z7 2 )- 1 - 2) = (det Vp) 2 + (det y y)2 - 2, 
since det (R Vip) = det V<p. In addition we note that £ H- (£ 2 + ^ — 2) is convex, hence 

( <detvrf2+ (5sW“ 2 ) 

is polyconvex and thus weakly lower semicontinuous. Concerning the other summands in the energy (1481) . we 
observe that ||sym(f7 — 1)|| is convex in U and both curvature contributions are quadratic. Hence, 

[ W{p,R)dV < liminf [ W(pk,Rk)dV (55) 

Jn k ^°° Jn 

and the proof is complete. 

Remark 15. For the general case regarding boundary conditions, i.e. P\ Td = Po we modify the energy functional 
to include a boundary contribution which realizes a weak coupling between microrotations and deformation 
gradient. We add 



||(I7-1) x vfdS, 


(56) 
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where again the cross product of a matrix with a vector is taken row wise. Then the NP W-Korn estimate (EH) 
holds with appropriate changes including the boundary contribution and the existence proof carries over. Note 
that 



||(t/ — 1) x v\\ 2 dS = 0 


{U-t)xv = 0 


a.e. 


(57) 


which implies 


(U-1)t 


(R T Vip — l).r 


0 for all tangential directions r. 


(58) 


The advantage of EH is that we do not presuppose a direct knowledge about the coupling of microrotations and 
deformations at the boundary r e> C dCl. Moreover, it is in general not possible to prescribe 

_ T 

(R X7ip — 1).t = 0 for all tangential directions r at Tjy, 
since the condition V</?.t = R.t cannot be met if^7<p is length changing onTu. 
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7 Appendix 

As mentioned in the notation section, the vector product x x y of two vectors x, y G R 3 has a close connection with skew-symmetric 
matrices. For v G R 3 fixed, we have 

( 0 -u 3 v 2 \ 

U3 0 —vi ] w V w G R 3 

-v 2 ui 0 / 

and for any skew-symmetric A G R 3x3 , there is v G R 3 such that Aw = v x w. A more geometrical interpretation of the vector 
product is given by its characterization “the vector product of a, b G R 3 is perpendicular to a and b both and its length is given by 
the area of the parallelogram spanned by a and 6”, which implies that the volume of a parallelepiped with sides a, b, c is given by 
(a x 6, c), as expressed by the following lemma: 

Lemma 16. Let a, 6 G R 3 . Then a x b is the unique element of R 3 which satisfies 

(a x 6, v) = vol(a, 6, v) for all v G R 3 (59) 


Proof. Instead of giving componentwise calculations let us recall the elementary geometric definition of the vector product: “The 
vector product of a, b G R 3 is perpendicular to a and b both and its length is given by the area of the parallelogram spanned by a 
and 6”, which implies that the volume of a parallelepiped with sides a, 6, c is given by {a x b,c). □ 

Lemma IT. Let u,v G R 3 and A G R 3x3 with adjugate Adj A. Then we have 

vol (Au, Av, Aw) = det A vol(u, v, w), 
vol (Au, Av, w) = vol(u, v, Adj Aw). 


Proof. Let u,v,w G R 3 and A G R 3x3 . According to the definition of vol (Definition [71) we have 

vo\(Au, Av, Aw) = det(Au\Av\Aw) = det j^4('u|u|u;)J 

= (det A) det(u|u|u;) = (det A) vol(u, v, w). 

Using this identity, the multilinearity of vol and the well known identity A Adj A = det A 1 we get 


vol (Au, Av, w) = 


det A 
1 


vol (Au, Av, det A w) 

det 

vol(Au, Au, A Adj Aw) = -vol(w, v, Adj Aw). 

det A 


Lemma 18. Let u,v,w G 1 


det A 

and AgR 3x3 . Then 

vol(Aw, v, w) + vol(u, Av, w) + vol(u, v, Aw) = (tr A) vol(w, v, w). 


□ 
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Proof. 


vol(Aix, v, w) + vol(ix, Av, w) + vol(ix, v, Aw) 
( a\\u\ + a\2U2 + 0 - 131 x 3 

= det I 021IX1 + 022 u 2 + O23IX3 


CL 31 U 1 + CL 32 U 2 + CL33U3 


Vl 

Wi\ 

/IX1 

“llfl + -212^2 + “ 13^3 

W 1 

V 2 

W 2 

+ det I 1X2 

“21^1 + 1222^2 + CL23V3 

W2 

V 3 

W 3 j 

\ix 3 

a 31U1 + II32V2 + 0332-3 

W 3 


III 

+ det I 1x2 

l w 3 


1 X 1 

1 x 2 

1>3 


0111/xi ■ 

0211/xi 

0311/xi 


- 0 , 121^2 ■ 

- d22W2 ■ 

- 03211 x 2 ■ 


- ai3U»3\ 

- d23W3 I 

- O33IXX3 / 


= dnU±V2W3 + di2U2V2W3 + 0131 x 31 x 211 x 3 + d2lU\V3W\ 

+ CL 22 U 2 V 3 W 1 + d 23 U 3 V 3 Wi + d 3 iU\V\W 2 + d 32 U 2 V\W 2 - 

— dnU\V 3 W 2 ~ d\ 2 U 2 V 3 W 2 ~ CI13U3V3W2 ~ d2\U\ViW3 


■ O33IX3IX1IIX2 


— d22U2VlW3 — O23O3IX1UX3 — CI31IX1IX2IIX1 — 0321X21X211X1 — O33IX3IX2UX1 
+ dnU 3 V\W 2 + d\2U3V2W2 + CL13U3V3W2 + d2\U\ViW3 

+ O22O1IX2UX3 + 0231X11X311X3 + O31IX2IX1UX1 + d 32 U 2 V 2 W\ + O33IX2IX3UX1 


— dnU2V\W3 — di2U2V2W3 — O13IX2IX3UX3 — O21IX3IX1IIX1 


— d 22 U 3 V 2 W\ — O23IX3IX3IIX1 — d 3 \U\ViW 2 ~ O32IX1IX2UX2 ~ O33IX1IX3IIX2 
+ di\U 2 V 3 Wi + di 2 U 2 V 3 W 2 + O13IX2IX3UX3 + O21IX3IX1UX1 

+ d 22 U 3 V\W 2 + d 23 U 3 V\W 3 + O31IX1IX2IIX1 + O32IX1IX2UX2 + O33IX1IX2UX3 


— dHU 3 V 2 Wl — dl 2 U 3 V 2 W 2 ~ O13IX3IX2IIX3 — O21IXIIX3IIXI 

— d22U\V3W2 — d23^\V3W3 — O31IX2IX1IIX1 — 0321X21X1 W2 ~ O33IX2IX1IXX3 

= 011 ( 1 x 11 x 211 x 3 — 1 x 11 x 311 x 2 + U 3 V 1 W 2 — U 2 V 1 W 3 + 1 x 21 x 311 x 1 — 1 x 31 x 211 x 1 ) 

+ 022 ( 1 x 2 ^ 3^1 — 1 x 21 x 11113 + 1 x 11 x 21113 — U 3 V 2 W 1 + U 3 V 1 W 2 — U 1 V 3 W 2 ) 

+ 033 (1x31x1 1x12 — 1x31x211x1 + 1x21131x11 — 1x11131x12 + U1V2W3 — U2V1W3) 
= (trA) det(ix |ix \w) 


□ 
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